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Abstract 


In  order  to  use  effectively  radiation  from  Josephson  junctions  Shitov  has  proposed  a  new 
design  to  couple  them  using  a  modified  mixed  series-parallel  array-biasing  scheme.  We  de¬ 
signed  two  different  layouts  of  samples  according  to  Shitov’s  proposal  and  performed  dc 
measurements  and  radiation  detection  experiments  at  the  frequency  of  80  GHz.  For  the  high 
temperature  range  we  observed  rather  large  radiation.  Furthermore,  we  studied  full  wave 
length  loop  antenna  arrays.  The  arrays  with  two  cells  have  more  complicated  dynamics,  simi¬ 
lar  to  discrete  rotobreathers  in  Josephson  ladders.  We  have  measured  detailed  current- voltage 
characteristics  of  single  and  double  loop  arrays  with  4x3,  4x4,  4x5,  4x6  and  4x7  junctions  per 
cell  in  the  broad  temperature  range  between  4.2  K  and  8  K.  In  order  to  decrease  the  space 
occupied  by  the  arrays  of  Josephson  junctions  studied  by  us  earlier,  we  replaced  bias  resistors 
by  Josephson  junctions.  We  designed  a  new  set  of  arrays  combining  both  the  junction  biasing 
and  resistor  biasing  schemes  on  the  same  chip.  The  comparison  of  two  alternative  biasing 
schemes  shows  that  the  junction  biasing  introduces  several  unwanted  complications  in  array 
characteristics.  Finally,  we  report  here  our  experimental  study  of  novel  qubit  devices.  The 
devices  have  been  designed  and  fabricated  using  Nb-A10a;-Nb  technology.  The  read-out  of 
qubit  state  has  been  performed  by  dc-SQUID.  We  observe  clear  signatures  of  coupling  be¬ 
tween  the  read-out  SQUID  and  a  qubit.  The  SQUID  response  has  a  hysteresis  related  to 
different  states  of  the  qubit. 


Chapter  1 

Series-parallel  arrays 


1  Layout  1 

Mixed  series-parallel  Josephson  junction  arrays  have  been  proposed  by  Jain  et  al  [1]  for 
high-frequency  oscillators.  Recently  Shitov  [2]  suggested  a  modification  of  that  scheme  using 
microwave  matching  of  dc  bias  leads.  A  sketch  of  such  an  array  consisting  of  only  5  junctions 
is  shown  in  Fig.l. 


Figure  1.1:  Series-parallel  Josephson  junction  array  scheme 

The  main  idea  of  this  layout  is  that  Josephson  junctions  (indicated  by  crosses  in  Fig.  1.1) 
are  arranged  in  parallel  for  DC  current  Ib  applied  to  the  array,  while  the  AC  voltages  add 
up  in  series.  The  interference  between  the  junctions  can  be  tuned  by  the  externally  applied 
magnetic  flux  through  the  cells.  At  half  frustration  the  amplitude  of  AC  oscillations  will  be 
the  largest  as  well  as  will  be  the  radiation  power.  The  DC  component  of  the  voltage  across 
the  dipole  antenna  attached  to  the  edge  junctions  of  the  array  should  be  equal  to  zero  for 
even  number  of  junctions  and  to  a  single  junction  voltage  for  odd  number  of  junctions.  The 
amplitude  of  AC  voltage  is  expected  to  increase  with  the  number  of  junctions.  The  array  has 
to  be  operated  in  a  resonant  mode  which  frequency  is  determined  by  the  cell  inductance  arid 
junction  capacitance. 
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2  Measurements  and  discussion:  Layout  1 

The  samples  according  to  our  layout  were  prepared  at  the  Hypres  foundry.  The  size 
of  each  junction  is  9 five?.  The  primary  cell  (containing  2  junctions  and  inductance)  size  is 
about  370  /xm^.  The  size  of  the  adjacent  (secondary)  cell  formed  by  the  leads  (vertical  lines 
in  Fig.  1.1)  and  lines  which  shorten  them  is  about  2980 ^m^.  We  have  designed  arrays  with 
10  junctions  and  with  20  junctions.  At  the  edges  of  the  array  we  put  capacitances  equal  to 
the  capacitance  of  a  junction.  The  arrays  were  designed  to  radiate  in  the  range  95-105  GHz. 
However  the  radiation  range  was  found  to  be  shifted  towards  lower  frequencies,  because  the 
critical  current  density,  dielectric  constant  and  other  parameters  came  out  slightly  different 
then  ones  which  were  taken  during  the  designing.  In  Fig.  1.2  we  present  the  layout  of  the 
array  of  20  junctions.  The  long  side  of  the  secondary  cells  is  about  a  quarter  wave  length  at 
the  center  frequency  (100  GHz).  This  means  that  the  loop  at  shorted  one  end  will  work  as 
open  at  this  frequency.  The  biasing  leads  of  the  array  end  with  chock  filters,  which  filter  noise 
in  the  frequency  range  of  80-110  GHz.  The  edges  of  the  array  are  attached  to  conductors  that 
form  a  dipole  antenna.  This  antenna  also  ends  with  chock  filters.  These  filters  are  needed 
to  prevent  the  leakage  of  the  radiation  outside  a  waveguide.  The  dipole  antenna  couples 
radiation  to  the  waveguide. 


2 


All  measurements  were  done  at  T  ~  7  K.  The  temperature  was  estimated  from  the  gap 
voltage  of  the  array,  which  was  about  1.3  mV.  To  change  the  frustration  through  the  array 
cell,  a  magnetic  field  perpendicular  to  the  array’s  plane  was  used.  As  usual,  first  we  did 
dc  measmements  of  the  current-voltage  characteristics  and  IdH)  pattern.  The  I-V  curve 
of  the  array  with  10  junctions  is  shown  in  Fig.  1.3(a)  by  dots.  This  curve  corresponds  to 
the  frustration  through  small  cells  of  f=0.4.  In  the  Fig.  1.3  the  dependence  of  the  critical 
current  versus  perpendicular  magnetic  field  is  presented.  This  curve  does  not  look  regular. 
We  suppose  that  it  is  due  to  the  large  secondary  cells,  which  are  about  eight  times  as  large 
as  the  primary  ones.  However,  it  still  shows  sufficiently  periodic  behavior.  In  the  presence  of 
an  external  frustration  in  larger  cells  we  may  expect  an  additional  supercurrent  which  brakes 
symmetry  of  the  system.  Moreover,  during  cooling  down  it  is  hard  to  avoid  trapping  random 
number  of  fluxons  in  larger  cells.  To  eliminate  this  problem  one  can  use  resistive  bridges 
instecid  of  superconducting  ones  (green  lines  in  Fig.  1.1),  but  in  this  case  a  considerable  part 
of  the  AC  current  would  dissipate  on  these  resistors. 


frustration 


Figure  1.3:  (a)  Current  voltage  characteristics  (dots)  and  measured  radiation  power  (line)  at 
80  GHz;  (b)  critical  current  pattern  versus  perpendicular  magnetic  field. 


Prom  the  I-V  curve  of  the  array  one  can  see  that  there  is  a  resonance  at  the  voltage  of 
about  180  fiV.  We  have  measured  radiation  from  this  array  using  our  room  temperature  su¬ 
perheterodyne  receiver.  By  black  line  in  Fig.  1.3(a)  we  show  the  measured  power  dependence 
on  voltage  across  the  array  for  the  presented  TV  curve.  On  this  curve,  there  is  a  very  narrow 
peak  at  the  same  voltage  as  the  resonance  branch  of  the  I-V  curve. 

The  weak  magnitude  of  the  detected  power  can  be  partially  explained  by  fault  design  of 
chock  filters.  All  chock  filters  at  the  end  of  dipole  antenna  are  also  part  of  this  antenna.  As 
a  part  of  the  antenna  they  should  be  placed  inside  the  waveguide.  The  waveguide  position  is 
shown  in  Fig.  1.2  by  blue  lines.  In  such  position  the  bonding  wires  would  be  most  probably 
working  as  RF  shorts  to  the  ground.  In  order  to  prevent  this  shorts,  we  had  to  make  the 
waveguide’s  walls  as  thin  as  it  is  shown  by  dashed  blue  lines.  This  complication  can  be  avoided 
by  moving  contact  pads  away  from  the  waveguide  walls. 

Thus,  in  these  experiments  with  Layout  1  the  radiation  was  measured  from  all  prepared 
samples.  The  detected  radiation  frequency  is  about  80  GHz.  We  note  several  design  improve¬ 
ments  to  be  made.  First,  we  would  have  to  avoid  additional  superconducting  loops.  The 
contact  pads  have  to  be  placed  away  from  the  chock-filters. 
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3  Layout  2 

In  the  Layout  2  we  improved  design  of  our  series-parallel  arrays.  First  improvement  in  this 
design  is  that  we  have  interrupted  the  secondary  superconducting  loops  by  resistor  bridges. 
The  bias  current  now  flows  through  the  bias  resistors.  Second,  the  capacitors  at  the  extremi¬ 
ties  of  the  array  have  been  substituted  by  Josephson  junctions.  Finally,  the  contact  pads  have 
been  moved  away  from  the  chock  filters,  so  that  the  bonding  wires  do  not  anymore  touch  the 
waveguide. 

4  Measurements  and  discussion:  Layout  2 

The  samples  according  to  our  micro-layout  were  prepared  at  the  Hypres  foundry.  The 
cell  size  of  all  prepared  arrays  is  about  350 /lim^.  The  length  of  the  inductance  leads  is  about 
300  jum.  We  have  designed  arrays  with  10  and  22  junctions.  The  size  of  junctions  was  either 
5x5  nm  2  or  3  X  3  /im  in  order  to  decrease  their  McCumber  parameter,  the  arrays  have  been 
measured  in  helium  vapor  at  temperatures  above  4.2  K.  The  actual  measurement  temperature 
was  estimated  from  the  gap  voltage.  Measurements  with  different  arrays  were  done  at  slightly 
different  temperatures.  The  parameters  of  the  arrays  and  experimental  data  are  summarized 
in  Table  1.1. 


Table  1.1:  Parameters  of  the  measured  samples 


Sample 

H) 

Number  of 
junctions 

Junction 
size,  /im^ 

Critical  current 
density,  A/cw? 

Measurement 
temperature,  K 

Measured  power 
level  (a.u.) 

C79#l 

22 

3x3 

910 

7.2 

1.0 

C79#2 

22 

5x5 

910 

7.5 

0.2 

C79#3 

10 

5x5 

910 

7.8 

0.0 

C92#l 

22 

3x3 

970 

7.8 

0.0 

C92#2 

22 

5x5 

970 

7.9 

57.0 

In  order  to  change  the  frustration  through  the  array  cell,  an  external  magnetic  field  per¬ 
pendicular  to  the  array’s  plane  was  used.  As  usual,  first  we  did  dc  measurements  of  the 
current- voltage  characteristics  and  /c(JT)  pattern.  From  dc  measurements  we  estimate  at 
which  frequencies  we  the  radiation  may  be  expected.  In  the  following  we  present  measure¬ 
ment  results  for  two  arrays.  Both  arrays  have  the  same  number  of  junctions  (22),  but  have 
different  junction  sizes. 

4.1  Sample  C79#l 

The  measurements  with  this  array  were  done  at  T=7.2  K.  The  dc  characteristics  of  the 
array  are  presented  in  Fig.1.4.  The  I-V  curve  is  shown  by  dots  in  Fig.  1.4(a).  This  curve 
corresponds  to  the  frustration  through  the  cells  of  f=0.23.  In  Fig.  1.4(b)  the  dependence  of 
the  critical  current  versus  magnetic  field  is  presented. 

Prom  the  I-V  curve  of  the  array  one  can  see  that  there  are  two  resonances  at  the  voltage 
of  about  87 /rV  and  260  |iV.  The  frequencies  of  both  resonances  (42  and  126  GHz)  are  out  of 
the  working  range  of  our  receiver  (78-120  GHz).  We  have  tried  to  measure  radiation  from  the 
array  at  frequencies  corresponding  to  the  harmonics  of  the  major  resonances.  At  the  frequency 
of  the  local  oscillator  (LO)  of  about  82  GHz  we  succeeded  to  observe  some  radiation.  The 
black  line  with  triangles  in  Fig.l.4(a)  corresponds  to  the  measured  radiation  power.  For  the 
LO  frequency  is  of  82  GHz,  the  radiation  peaks  are  located  near  to  the  voltages  of  the  lower 
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Figure  1.4:  dc  characteristics  of  the  array  C79#l:  (a)  I-V  curve  (dots)  and  measured  power 
vs.  voltage  (triangles),  (b)  critical  current  vs.  frustration  pattern 


resonance,  which  corresponds  to  42  GHz.  We  believe  that  the  detected  radiation  is  the  second 
harmonic  of  the  lower  resonance. 

In  Fig.  1.5 (a)  we  present  the  dependence  on  the  frustration  for  the  critical  current  (dots) 
and  maximum  measured  radiation  power  (triangles).  The  power  vs.  frustration  curve  looks 
quite  regular.  Prom  this  curve  we  deduce  that  the  array  radiates  maximum  power  at  82  GHz 
for  frustration  around  0.25  and  0.75.  Fig.l.5(b)  presents  I-V  curves  and  measured  radiation 
power  (inset)  at  frustration  0.23  (black  symbols)  and  0.78  (red  line).  The  data  corresponding 
to  the  frustration  of  0.78  are  asymmetric.  At  the  positive  voltages  two  I-V  curves  look  the 
same,  power  vs.  voltage  curves  are  also  very  similar.  I-V  curves  at  the  negative  voltages  are 
slightly  different,  consequently  the  measured  power  is  also  different. 


Figure  1.5:  (a)  Frustration  dependence  of  the  critical  current  (dots)  and  measured  power 
(triangles)  of  the  sample  C79#l;  I-V  curves  and  measured  power  (inset)  at  frustration  0.23 
and  0.78 
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4.2  Sample  C92#2 

The  measurements  with  this  array  were  done  at  a  little  higher  temperatures  (7.9  K). 
The  measured  I-V  curve  and  Ic{H)  pattern  are  presented  in  Fig.  1.6.  On  the  I-V  curve  (dots 
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Figure  1.6:  dc  characteristics  of  the  array  C92#2:  (a)  I-V  curve  (dots)  and  measured  power 
vs.  voltage  (triangles),  (b)  critical  current  vs.  frustration  pattern 


in  Fig.l.6(a))  of  the  array  there  is  resonance  at  85  fiY.  The  higher  resonance  (which  was 
observed  at  lower  temperatures  at  voltage  of  about  165  //V)  is  suppressed  at  this  temperature 
due  to  high  dissipation.  However,  we  were  able  to  measure  relatively  strong  radiation  from 
the  array  close  to  the  voltage  corresponding  to  the  higher  resonance.  The  LO  frequency  was 
set  to  80  GHz.  The  triangles  in  Fig.  1.6(a)  show  the  measured  power.  This  curve  has  two 
sharp  peaks  around  the  voltage  of  165  /zV.  In  the  Fig.  1.6(b)  the  dependence  of  the  critical 
current  versus  magnetic  field  is  presented. 

In  Fig.  1.7  we  present  the  critical  current  (dots)  and  the  radiation  power  (triangles)  depen¬ 
dence  on  the  frustration.  One  can  see  that  the  radiation  is  at  maximum  at  nearly  half-integer 
frustration.  Small  detuning  of  the  frustration  sharply  reduces  radiation.  The  maximum  mea¬ 
sured  power  was  about  sixty  times  large  than  that  of  the  sample  079^1.  This  is  partially 
due  to  the  fact  that  the  main  harmonics  have  higher  radiation  power  and  at  least  partially 
due  to  the  larger  size  of  the  junctions  of  the  array  C92#2  and,  consequently,  larger  dc  power 
supplied.  We  suppose  also  that  better  coupling  of  the  array  to  the  waveguide  was  achieved 
in  latter  case. 

Thus,  from  sample  C79#l  we  have  detected  radiation  at  frequency  of  82  GHz,  which 
corresponds  to  the  second  harmonic  of  the  lower  resonance.  The  largest  power  was  detected 
from  the  sample  C92#2.  The  frustration  dependence  of  the  measured  power  is  found  to  be 
quite  regular  for  all  arrays. 
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Figure  1.7:  Frustration  dependence  of  the  critical  current  (dots)  and  measured  power  (trian¬ 
gles)  of  the  sample  C92#2 
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Chapter  2 

Full-wavelength  loop  antenna  arrays 


1  Loop  antenna  array 

A  Josephson  junction  row  array,  where  the  perimeter  of  each  rectangular  cell  in  the 
row  of  about  one  wavelength,  may  provide  a  simple  way  to  achieve  a  phased  array  sub-mm 
radiation  source.  The  beam  formed  above  the  array  may  be  directed  by  imposing  a  uniform 
external  magnetic  field  over  the  whole  array  and  thereby  inducing  a  constant  phase  shift  in 
the  emitted  radiation  from  cell  to  cell.  The  main  concern  with  this  type  of  array  is  that 
due  to  fabrication  defects,  the  spread  in  the  DC  voltages  (and  consequent  frequencies)  of  the 
contributing  junction  oscillators  may  overwhelm  the  contribution  of  the  oscillating  antenna 
current  that  tends  to  keep  the  junction  oscillators  in  synch.  Preliminary  calculations  by  S. 
Yukon  and  F.  Lin  indicate  that  when  the  junctions  are  matched  to  the  antenna  impedance, 
these  currents  can  be  significant  and  may  allow  such  a  system  to  work  [1].  In  Fig.  2.1  we 
present  electric  scheme  of  single  (a)  and  double  (b)  cell  arrays  with  4x5  Josephson  junctions 
per  cell. 


Figure  2.1:  Equivalent  scheme  of  a  loop  antenna  array  with  one  (a)  and  two  (b)  cells 


2  Measurements  and  discussion 

According  to  our  layout  and  specification,  the  samples  having  the  critical  current  density 
Jc  ~  1000 A/cm^  were  prepared  at  the  foundry  HYPRES  [2].  The  samples  are  single  cell  or 
double  cell  arrays.  The  arrays  are  with  4x3,  4x4,  4x5,  4x6  and  4x7  around  a  cell.  The  size 
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of  each  junction  is  about  9  Fig.  2.2  shows  a  micro  photograph  of  the  sample  #34- 

1x4x5.  Table  2.1  presents  all  specific  parameters  of  the  prepared  samples.  The  inductance 
was  estimated  using  numerical  results  by  Ketchen  [3]. 


Figure  2.2:  Micro  photograph  of  the  sample  #34-1x4x5,  loop  antenna  array  with  4x5  junctions 


Table  2.1:  Estimated  parameters  of  the  samples 


Number  of 
junctions 

Cell  size, 
^ivc? 

Total 

capacitance,  fF 

Cell 

inductance,  pH 

Resonance 
frequency,  GHz 

4x3 

4356 

30 

147.6 

75.6 

4x4 

7056 

22.5 

194.2 

76.1 

4x5 

10404 

18 

249.9 

75.0 

4x6 

14400 

15 

301.6 

74.8 

4x7 

19044 

12.9 

355.5 

74.3 

The  layout  of  a  single  cell  array  with  4x5  junctions  is  presented  in  Fig.  2.3.  The  ex¬ 
ternal  bias  current  (/bias)  is  injected  at  each  node  of  the  array  via  67  O  resistor  (along  the 
vertical  direction  in  Fig.  2.3)  and  voltage  (F)  is  measured  across  the  vertical  junctions.  The 
description  of  the  used  measurement  setup  can  be  found  in  one  of  our  previous  reports  [4]. 


fcottom  Nbo 
top  NbD 
resistor  (Mo) 


Figure  2.3:  Layout  of  a  cell  with  4x5  junctions 

We  have  measured  current  voltage  (I-V)  characteristics  of  the  arrays  at  the  temperature 
range  4.2-8  K.  Due  to  their  large  inductance  the  arrays  are  insensitive  to  external  magnetic 
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field  perpendicular  to  the  array  plane. 

In  Fig.  2.4  we  present  several  I-V  curves  of  single  loop  arrays  taken  with  storage  oscil¬ 
loscope  at  different  temperatures.  I-V  characteristics  of  double  cell  arrays  taken  at  different 
temperatinres  are  presented  in  Fig.  2.5. 


(a)  4x4  junctions,  T=6.9  K,  chip  #2  (b)  4x5  junctions,  T=5.9  K,  chip  #1 


(c)  4x6  junctions,  1=7. Z  K,  chip  #1  (d)  4x6  junctions,  T=4.2  K,  chip  #2 


Figure  2.4:  I-V  characteristics  of  single  cell  arrays  with  4x4  (a),  4x5  (b),  4x6  (c)  and  4x7 
(d)junctions  around  a  cell 


Fig.  2.4(a)  presents  I-V  characteristics  of  a  cell  with  4x4  junction  measured  at  T=6.9 
K.  PYom  this  figure  we  note  that  all  of  four  vertical  junctions  switch  to  resistive  state  at 
nearly  the  same  bias  current.  However,  in  spite  of  their  small  spread  we  could  measure  all 
the  resistive  branches.  To  measure  n-th  branch  we  increase  (decrease)  bias  current  until  the 
system  switches  to  this  state  and  decrease  (increase)  it  until  the  system  switches  to  another 
branch.  On  this  I-V  curve  we  find  no  no  resonance  branches.  We  did  not  observe  any 
resonance  in  the  whole  temperature  range  that  we  measured. 
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(a)-(b)  4x5  junctions,  T=4.5  K,  chip  #2 


(c)  4x7  junctions,  T=4.5  K,  chip  #2 


(d)  4x7  junctions,  T=4.2  K,  chip  #1 


■ 

■ 

m 

■ 

■ 

1 

■ 

1 

K 

■ 

■ 

M 

■ 

■ 

1 

U 

■ 

■ 

■ 

■ 

M 

■ 

■ 

1 

■ 

■ 

■ 

i 

Hi 

■ 

1 

■ 

■ 

■ 

p 

■ 

1 

Mi 

■ 

■ 

r 

i 

'v, 

I 

■ 

■ 

■ 

1 

1 

■ 

r 

■ 

■ppr  A' 

ii 

Ml 

m 

■ 

■ 

i 

iB 

B 

Figure  2.5:  I-V  characteristics  of  double  cell  arrays  with  4x5  (a)-(b)  arid  4x7  (c)-(d)  per  cell 


Fig.  2.4(b)  shows  I-V  curve  of  a  cell  with  4x5  junctions  at  T=5.6  K.  This  sample  shows 
some  resonance  behavior.  The  first  resonance  step  (indicated  by  red  arrow  in  Fig.  2.4(b)) 
a,ppears  at  frequency  about  710  GHz,  which  is  much  higher  then  the  expected  firequency  (~75 
GHz).  The  same  steps  can  be  seen  at  the  bottom  part  of  all  branches.  One  can  also  see  that 
the  critical  currents  of  difierent  branches  are  different,  so  that  the  sample  is  not  uniform  and 
it  is  hard  to  interpret  its  behavior. 

Fig.  2.4(c)  presents  I-V  curve  of  a  single  cell  array  with  4x6  junctions.  The  measurement 
were  done  at  temperature  7.3  K.  This  sample  shows  nearly  the  same  critical  current  of  different 
branches,  n-th  branch  corresponds  to  the  state  with  n  vertical  junctions  are  in  resistive  state. 
We  also  did  not  observe  any  resonance  with  this  sample. 

The  I-V  curve  presented  in  Fig.  2.5(d)  corresponds  to  2  cell  array  with  4x7  junctions  per 
cell.  The  first  6  branches  are  similar  to  that  of  a  single  cell  array.  However,  the  7-th  branch 
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itself  splits  to  other  7  branches.  In  order  to  measure  them  we  increased  bias  current  until 
some  value  (we  denote  it  as  /stop)  ^d  than  decrease  it  until  zero.  Depending  on  /stop  different 
branches  appear  on  the  I-V  curve.  We  believe  that  they  correspond  to  different  configurations 
of  several  horizontal  and  vertical  junctions  in  resistive  state.  Such  a  states  are  called  discrete 
rotobreathers  and  have  been  observed  in  Josephson  leidder  arrays  [5,  6). 

We  have  measured  more  than  15  different  arrays.  All  arrays  with  single  cell  show  similar 
behavior  like  in  Fig.  2.4(a)-(c).  This  corresponds  to  one  by  one  switching  to  resistive  state 
of  vertical  junctions.  The  arrays  with  two  cells  have  more  complicated  dynamics.  They  show 
different  branches  not  only  at  the  bottom  part  of  the  last  gap,  but  also  at  the  bottom  parts 
of  other  gap  branches. 

3  Conclusion 

Measurements  with  full  wave  length  loop  antenna  arrays  have  been  performed.  We  did 
not  observe  any  systematic  resonances.  All  of  the  single  cell  arrays  have  shown  more  or 
less  uniform  critical  current  of  different  branches  on  the  I-V  curve.  The  arrays  with  2  cells 
has  shown  different  breather  states.  In  order  to  investigate  exactly  which  junctions  are  in 
the  resistive  state  at  the  given  voltage  one  may  make  additional  probes  to  measure  voltages 
across  each  junction.  An  easier  possible  method  is  using  our  low  temperature  scanning  laser 
microscope  [6]. 
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Chapter  3 

Replacement  of  resistors  by 
Josephson  junctions  in  biasing  leads 

1  Introduction 

The  conventional  biasing  scheme  of  Josephson  junctions  arrays  consist  of  resistor  network 
(Fig.3.1a).  In  order  to  have  a  homogenous  bias  current  over  the  whole  array,  the  resistances 
in  the  biasing  have  to  be  large  in  comparison  with  junction  sub-gap  resistance.  Such  resistors 
extend  over  several  cell  lengths  from  the  array  and  may  decrease  the  radiation  efficiency. 


Figure  3.1:  Schematic  top  view  of  the  arrays  with  resistor  biasing  (a)  and  jimction  biasing 
(b)  networks.  Red  junctions  belong  to  active  oscillator  row,  while  green  junctions  form  the 
bias  network.  Dashed  lines  show  the  possible  flux  flow  direction  through  the  arrays. 

Following  recent  proposal  by  Stanford  Yukon  [1],  we  have  replaced  the  bias  resistors  by 
Josephson  junctions  (Fig.  3.1b).  According  to  theoretical  calculations,  the  radiation  efiiciency 
is  supposed  to  be  around  3-4%  larger  than  that  of  arrays  with  resistive  biasing.  Because  the 
bias  junctions  are  supposed  to  switch  to  the  resistive  state  at  a  cmrent  equal  or  less  than  the 
switching  current  of  the  vertical  junctions  in  the  active  row,  one  iriay  neglect  the  Meissner 
state  in  the  biasing  network. 

In  case  of  resistive  bias,  the  resistor  connected  to  the  corner  vertical  junction  is  twice 
larger  than  resistors  connected  to  two  vertical  junction.  Following  this  idea,  the  bias  junc¬ 
tion  connected  to  only  one  vertical  junction  should  be  twice  smaller  than  the  bias  junction 
supplying  current  to  two  vertical  junctions  of  the  active  array. 


15 


2  Array  layout 

The  design  of  the  arrays  schematically  shown  in  Fig.  3.1  was  based  on  our  previously  studied 
arrays  [2].  The  active  row  junctions  have  an  area  of  9/xm^  and  the  cell  size  is  150  )um^.  The 
junction  area  of  9  fj.va?  is  close  to  the  smallest  permitted  by  the  Hypres  design  rules  and  we 
decided  to  keep  it  unchanged  from  the  previously  made  and  tested  design.  The  bias  network 
have  junctions  with  an  area  of  18  /xm^  (edge  junction  of  9  ^m^).  The  total  area  (and  thus,  the 
field-free  critical  current)  of  each  of  the  bias  rows  is  equal  to  that  of  the  active  row.  Thus,  at 
zero  frustration  and  absence  of  self-field  effects  all  rows  are  supposed  to  switch  to  the  resistive 
state  at  the  same  bias  current. 

In  order  not  to  affect  the  resonance  properties  of  the  active  row  the  bias  cell  size  was 
chosen  much  larger  than  the  size  of  active  row  cells.  We  have  designed  arrays  with  two 
different  cell  size,  3000 /xm^  (sample  #1)  and  700  ^m^  (sample  #2).  As  a  reference  array  we 
have  put  an  array  with  resistor-based  biasing  (reference  sample).  The  resistors  are  about 
25  fl  (edge  resistor  ~50f2).  All  designed  arrays  have  12  cells  in  the  active  row.  All  horizontal 
junction  are  embedded  into  a  stripline  formed  between  the  array  electrodes  and  the  ground 
plane  below  them. 

3  Experimental  results 

According  to  our  design,  the  samples  were  fabricated  at  Hypres  foundry.  The  critical 
current  density  of  the  junctions  is  about  1000  A/cm^  at  4.2  K.  All  measurements  were  done 
at  7.9  K.  We  have  systematically  measured  current-voltage  characteristics  and  perpendicular 
magnetic  field  Ic(H)  pattern  of  the  critical  current  and  the  resonance  voltage. 


Figure  3.2:  Magnetic  field  pattern  of  the  critical  current  (a)  and  current-voltage  characteristics 
at  f=0.3  (b)  of  the  array  with  resistor  based  bias  network. 

In  Figure  3.2a  we  present  Ic(H)  pattern  of  the  reference  array  with  biasing  via  resistors.  As 
expected,  we  find  regular  oscillations  of  the  critical  current.  Its  current  voltage  characteristics 
at  frustration  f=0.3  shows  well  pronounced  resonance  at  voltage  of  about  170 /xV. 

The  dependence  of  the  critical  current  of  the  active  row  of  sample  #1  is  shown  Fig.  3.3a. 
For  comparison,  we  present  the  Ic(H)  pattern  of  the  bias  row  in  Fig.  3.3b.  In  spite  of  the 
regular  oscillation  of  the  bias  rows’  critical  current,  the  critical  current  of  the  active  row  is 
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Figure  3.3:  Magnetic  field  patterns  of  the  critical  current  of  the  active  row  (a)  and  of  the  bias 
rows  (b)  of  the  array  with  junction-based  biasing.  Biasing  cell  size  is  3000 /im^. 


rather  irregular.  Note,  that  in  order  to  exclude  possible  effects  of  residual  fields  the  reference 
array  was  measured  in  the  same  run  as  samples  with  junction-based  biasing  network.  Second 
peculiar  point  is  that  the  active  row  critical  current  is  smaller  than  that  of  biasing  rows.  The 
ratio  of  the  oscillation  periods  of  the  bias  rows  and  active  row  is  about  three  times  smaller 
than  the  ratio  of  the  corresponding  cell  sizes.  The  critical  currents  of  the  bias  row  also  differ 
from  one  another. 

The  I-V  curve  of  sample  #1  is  presented  in  Fig.  3.4  at  frustration  f=0.4.  The  I-V  curve 
of  the  active  row  is  similar  to  that  of  the  array  with  resistive  biasing.  The  I-V  curve  of  the 
second  biasing  row  (red  symbols)  shows  a  small  peculiarity  by  deviating  from  zero  voltage 
state  at  bias  current  of  about  90 /xA.  This  biasing  row  remains  in  the  resistive  state  up  to  the 
current  Ibias=110AiA  and  then  returns  back  to  the  superconducting  state. 


Figure  3.4:  Current  voltage  characteristics  of  the  array  with  junction  based  biasing.  Biasing 
cell  size  is  3000 /xm^. 

This  peculiar  reentrant  behavior  has  been  previously  observed  by  us  for  symmetry-broken 
states  in  Josephson  junction  arrays.  It’s  spatial  patterns  have  been  visualized  by  using  low 
temperatm-e  scanning  laser  microscope  [4,  5].  Possible  explanation  for  the  array  dynamics  in 
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the  measurements  presented  here  is  illustrated  in  Fig.  3.1.  The  junctions  marked  by  circles 
are  in  resistive  state  and  the  dashed  line  shows  how  the  magnetic  flux  travels  across  the  array. 
Thus,  symmetry-broken  states  occur  when  one  or  more  horizontal  junctions  (placed  in  the 
direction  perpendicular  to  the  bias  current)  go  resistive. 


Figme  3.5;  Magnetic  field  patterns  of  the  critical  current  of  the  active  row  (a)  and  of  the  bias 
rows  (b)  of  the  array  with  junction-based  biasing  scheme.  Biasing  cell  size  is  700 /im^. 

The  dependence  of  the  critical  current  of  the  active  row  of  sample  #2  is  shown  Fig.  3.5a 
and  the  Ic(H)  pattern  of  the  bias  rows  in  Fig.  3.5b.  These  patterns  are  similar  to  that  of 
sample  #1.  The  ratio  of  the  areas  of  the  cells  in  the  biasing  and  active  rows  is  about  4. 
However,  the  oscillation  period  of  the  critical  current  of  the  bias  rows  is  same  as  that  in  the 
active  row. 

The  I-V  curves  the  array  #2  (Fig.  4.4,  f=0.7)  are  similar  to  the  I-V  curve  of  the  sample 
The  peculiar  reentrant  resistive  state  appears  now  for  both  bias  rows  (red  and  blue 
symbols),  but  at  different  bias  currents. 


Voltage  (nV) 


Figure  3.6:  Current-voltage  characteristics  of  the  array  with  the  junction-based  biasing 
scheme.  Biasing  cell  size  is  3000 /tm^. 
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Finally,  by  biasing  arrays  at  constant  bias  current,  we  have  measured  the  magnetic  field 
patterns  of  the  voltage  at  the  resonant  step.  In  Fig.3.7  we  present  F(H)  patterns  at  Ibias  = 
120 /xA  for  the  all  three  arrays  are  presented.  Similar  to  critical  current  patterns,  the  K(H) 
patterns  of  the  arrays  with  junction  biasing  are  irregular. 


Frustration 


— -S=700  Mm2 
—  S=3000  um2[ 


0.2  0.4  0.6  0.8  1.0  1.: 

Frustration 


Figure  3.7:  Volt-field  pattern  of  the  reference  array  (a)  and  of  the  arrays  with  junction-based 
biasing  scheme  (b). 


4  Conclusion 

We  have  experimentally  studied  the  array  of  Josephson  junctions  with  biasing  via  a  row 
of  Josephson  junctions.  We  have  observed  that,  in  the  operational  range  of  firustration,  the 
active  oscillator  row  switches  to  the  resistive  state  at  lower  current,  at  which  biasing  junctions 
are  still  in  the  superconducting  state.  We  also  find  strong  unwanted  dynamic  interaction 
between  the  oscillator  and  biasing  rows  of  junctions.  This  dynamic  interaction  most  probably 
corresponds  to  symmetry  broken  states.  ■- 

The  area  (and  thus  the  critical  ciurrent)  of  biasing  junctions  in  this  layout  was  limited 
by  Hypres  design  rules.  In  principle,  by  modifying  the  design  of  our  standard  oscillator  to 
larger  junctions  we  should  be  able  to  make  biasing  junctions  switching  before  the  active  row. 
This  would  require  additional  test  run  in  order  to  vary  the  modified  cell  sizes  and  match  the 
radiation  frequency  range  of  the  receiver  at  80-120  GHz.  Nevertheless,  we  suppose  that  the 
dynamic  interaction  between  the  oscillator  and  the  biasing  rows  will  remain  a  serious  problem 
in  such  a  design. 


19 


Bibliography 

[1]  Stanford  Yukon,  private  communication  (2002) 

[2]  P.  Caputo  and  A.V.  Ustinov,  Iterim  report  on  ’’Analjreis  of  Triangular  Arrays  of  Joseph- 
son  Tunnel  Junctions”,  contract  F61775-01-W-E041  (1998) 

[3]  http://www.liyiires.com 

[4]  D.  Abraimov,  P.  Caputo,  G.  Filatrella,  M.  V.  Fistul,  G.  Y.  Logvenov,  and  A.  V.  Ustinov, 
Phys.  Rev.  Lett.  83,  5354  (1999). 

[5]  P.  Binder,  D.  Abraimov,  A.  V.  Ustinov,  S.  Flach  find  S.  Zolotaryuk,  Phys.Rev.Lett.  84, 
pp.  745-748,  (2000) 


20 


Chapter  4 

Triangular  prism  qubits 


1  Introduction 

Several  novel  qubit  devices  have  been  proposed  by  Stanford  Yukon  [1,  2].  The  most 
intriguing  of  these  is  a  discretized  version  of  the  triangular  long  Josephson  junction  prism 
[3-7].  Currently,  there  is  no  technical  method  to  make  such  a  real  3D  system.  There  has 
been  proposed  [8]  a  flattened  2D  version  of  the  qubit,  which  is  shown  in  Fig.  4.1.  Potential 
energy  proflle  of  this  qubit  is  periodic  and  each  period  has  two  minima.  These  minima  of 


Figure  4.1:  2D  version  of  Yukon’s  qubit  and  a  read-out  SQUID. 


the  potential  profile  of  this  system  are  associated  with  the  circulating  current  in  the  left 
(ii  =  imin,  in  =  0)  and  right  loops  (i/,  =  0,  iji  =  Imin)-  The  flux  and  -f-  determine 
the  barrier  height.  With  change  of  barrier  height  its  the  width  of  the  wells  and  of  the  barrier 
and  /min  change.  By  changing  the  differential  flux  between  the  flux  through  the 

left  and  right  loops  one  can  make  the  potential  profile  asymmetric.  The  qubit  state  can  be 
measured  with  help  of  read-out  SQUID  placed  next  to  it. 

The  standard  qubit  based  on  the  single-loop  SQUID  has  a  double  well  potential  at  half 
frustration  and  relies  on  the  external  field.  The  important  feature  of  the  Yukon’s  qubit  is  that 
it  should  be  less  sensitive  to  fluctuations  of  the  external  magnetic  field.  Another  interesting 
property  of  this  qubit  is  that  the  barrier  height  can  be  changed  without  changing  the  symmetry 
of  the  potential  profile.  The  main  difference  between  multi-qubit  circuit  based  on  the  Yukon’s 
qubits  and  single-loop  qubits  is  that  former  ones  can  be  coupled  by  Cirac-Zoller  tsrpe  of  bus. 

[9] 
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2  Design 

In  order  to  characterize  the  potential  energy  profile  in  the  classical  regime,  we  have  made 
design  of  the  circuit  shown  schematically  in  Fig.  4.1.  There  are  several  constrictions  for  the 
parameters  of  the  qubit  and  dc-SQUID.  These  constrictions  are  set  by  the  design  rules  of  the 
foundry  [10]. 

In  the  theoretical  calculations  the  qubit  cell  inductance  is  neglected.  However,  technolog¬ 
ical  processes  put  the  lower  limit  on  the  size  of  cell.  In  order  to  minimize  the  cell  inductance 
of  the  qubit  we  made  the  cell  size  2x2  /rm.  The  roughly  estimated  inductance  of  the  cell  is 
about  L  =  3  pH.  The  junction  size  is  about  3x3  jam. 

For  making  standard  current-biased  SQUID  read-out,  the  junctions  of  the  SQUID  should 
be  overdamped  and  its  fii,  parameter  should  be  of  the  order  of  one.  In  order  to  have  an 
overdamped  case,  junctions  were  shunted  by  a  resistance  of  5  f2.  The  cell  size  of  the  SQUID 
was  made  8x8 /am.  The  layout  is  shown  in  Fig.  4.2. 


Figure  4.2;  Design  layout 


A  magnetic  flux  through  the  SQUID  cell  is  controlled  by  a  current  through  the  control 
line  (CL)  marked  as  Isqcl.  The  qubit  state  is  controlled  by  currents  through  CLs  marked 
as  II  and  Ir.  Unfortunately,  the  circuits  designed  with  a  middle  CL  (for  applying  differential 
flux)  came  defected  from  the  foundry  due  to  a  limitation  given  by  design  rules.  Samples  either 
with  or  without  ground  plane  were  designed  on  the  chip. 

3  Experiment 

Samples  according  to  our  micro  design  were  prepared  at  Hypres  foundry.  The  critical 
current  density  was  about  1000  A/cm^  at  4.2  K.  In  spite  that  the  junctions  were  shunted, 
the  read-out  SQUID  was  slightly  in  a  hysteretic  regime.  In  order  to  make  it  sufficiently 
overdamped  we  have  increased  temperature  up  to  5K.  As  expected,  the  samples  with  the 
ground  plane  were  almost  imsensitive  to  the  fields  produced  by  CLs  for  currents  up  to  50  mA. 

First,  we  have  characterized  the  read-out  SQUID.  In  Fig.  4.3  we  present  the  measured 
current-voltage  characteristics  (a)  and  critical  current  versus  frustration  produced  by  Isqcl 
(b).  Prom  this  curves  we  have  estimated  the  SQUID’s  parameters  ss  fit,  =  2.66  and  0c  —  0.25. 

As  a  working  point  of  the  read-out  SQUID  we  have  chosen  the  bias  current  and  frustration 
with  largest  dV/dl  and  dlc/dfsq  response  amplitudes.  This  point  is  marked  by  a  solid  dot  on 
I-V  curve  (1=  164/iA)  and  Ic(H)  pattern  (fsg  =  0.75).  We  have  fixed  the  current  Ir  through 
the  right  CL  while  sweeping  the  current  of  II  the  left  CL  and  measuring  the  voltage  of  the 
SQUID.  In  Fig.  4.4  we  present  voltage-field  curves  corresponding  to  increase  of  frustration 
(black  line)  and  decrease  of  frustration  (red  line)  for  Ir  =  153^A.  The  inset  presents  the 
difference  between  this  curves.  The  differential  flux  between  the  left  and  right  loops  vanishes 
($i  =  $jj)  at  II  =  153/iA  (/SQ  =  0.2). 
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Figure  4.3:  Current-voltage  characteristics  (a)  and  critical  current  pattern  vs.  frustration  (b) 
of  the  read-out  SQUID 

l^(mA) 

-0.70  -0.35  0.00  0.35  0.70  1.05  1.40 


Figure  4.4:  Read-out  SQUID  response  to  increase  (black  line)  and  decrease  (red  line)  of 
current  II-  Inset:  the  difference  between  the  two  curves. 

The  hysteresis  shown  in  Fig.  4.4  starts  at  fsQ  ci-0.86.  At  fsQ  ^0.3  (red  arrow)  it 
abruptly  increases  and  again  decreases  at  fsQ  iil.O  (black  arrow).  The  hysteresis  vanishes 
at  fsQ  ~1.98.  The  effective  change  of  the  frustration  through  SQUID  for  larger  hysteresis  is 
about  Af~15.3  xlO“®  and  for  smaller  one  it  is  Afe8.7  xlO“^.  We  suppose  that  large  hystere¬ 
sis  occurs  due  to  a  switch  of  the  system  between  the  left  (Fig.  4.5a)  and  the  right  (Fig.  4.5b) 
minima. 

The  smaller  hysteresis  occurs  due  to  change  of  sign  of  a  mesh  current  in  one  of  the  qubit 
cells.  For  example,  between  states  shown  in  Fig.  4.5b  and  Fig.  4.5c  the  mesh  current  in 
changes  its  direction.  In  order  to  better  understand  the  double-hysteretic  behavior  we  plan 
to  perform  another  experiment  with  CLs  which  change  the  flux  differentially  and/or  change 
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Figure  4.5:  For  explanation  of  two  hysteresis:  a)  sjretem  is  in  the  left  minimum;  b)  system  is 
in  the  right  minimum;  c)  system  is  in  right  minimum,  but  with  a  negative  mesh  current. 

the  total  flux  through  left  and  right  loop,  leaving  the  flux  in  the  middle  loop  unchanged. 

4  Conclusion 

We  have  performed  first  experiments  with  novel  Josephson  jimction  qubits.  We  have 
observed  that  the  read-out  response  has  a  double-hysteretic  behavior.  Next  step  will  be  to 
investigate  this  hysteresis  in  dependence  on  the  second  control  line  current  and  to  find  out 
the  magnetic  field  periods  of  the  qubit  by  its  dc  characterization.  We  note  several  design 
improvements  to  be  made.  First,  we  have  to  improve  the  coupling  of  our  read-out  system. 
Next,  a  differential  flux  control  line  has  to  be  added.  It  planned  to  develop  an  improved 
version  of  this  layout  in  the  near  future. 
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